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Abstract 

The purpose of this article is to extend the lightcone worldsheet lattice description of 
^^ ^ string theory to include the Neveu-Schwarz model. We model each component of the 

,— I- ■ fermionic worldsheet field by a critical Ising model. We show that a simple choice of 

CN ■ boundary conditions for the Ising variables leads to the half integer modes required by 

the model. We identify the G-parity operation within the Ising model and formulate 
the procedure for projecting onto the even G-parity sector. We construct the lattice 
^^ . version of the three open string vertex, with the necessary operator insertion at the 

^^ \ interaction point. We sketch a formalism for summing planar open string multi-loop 

amplitudes, and we discuss prospects for numerically summing them. If successful, 
the methods described here could provide an alternative to lattice gauge theory for 
KA , computations in large N QCD. 
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1 Introduction 

Nonabelian gauge theory in 4 spacetime dimensions can be regarded as the zero slope hmit 
{a' — )■ 0) of open string theory [1] in a higher dimension (e.g. D = 10), provided that 
the open strings are required to end on D3-branes [2]. For the conformally invariant A/" = 4 
supersymmetric gauge theory, the AdS/CFT correspondence asserts that, after reinterpreting 
the planar open string multiloop diagrams as closed string trees, the a' — )■ limit remains 
a string theory on the curved background AdS5xS^ [3]. In this sense A/" = 4 isn't just the 
limit of string theory: it is a (closed) string theory. Moreover, the conformal invariance of 
A/" = 4 means that the gauge coupling is a true parameter, and this closed string theory can 
be analyzed semiclassically in the limit of large 't Hooft coupling Ng'^. 

For QCD, however, asymptotic freedom precludes such a parametric semi-classical closed 
string limit. Even though 't Hooft's large A^ limit [4] of QCD should still have an inter- 
pretation in terms of some kind of closed string background, the absence of a semiclassical 
approximation to find and analyze that background casts doubt on the practical value of 
a closed string interpretation in actual calculations. We entertain here the possibility that, 
in this circumstance, it may be more profitable to forego the equivalence of QCD to some, 
as yet to be discovered, closed string theory and instead to exploit the known connection 
of large N QCD to open string theory, keeping a' > as a regulator until the end of the 
analysis. The hope is that some nonperturbative aspects of QCD, such as quark confinement 
and the meson spectrum in 't Hooft's large A^ limit, will be more tractable treated as an 
open string theory with a' > than as a quantum field theory, its a' — > limit. 

Several years ago my collaborators and I showed how to represent each gauge theory 
planar diagram as a lightcone open-string worldsheet path integral [5] . This work established 
a direct formal connection of large A^ QCD to an open string theory with a' = 0. The 
problem with keeping a' = 0, however, is that the UV divergences of quantum field theory 
are not properly regulated, so that various counterterms must be bought in, order by order 
in the loop expansion, to cancel UV induced artifacts that violate Lorentz invariance [6]. 
Keeping a' > certainly mitigates these problems, although it remains to be seen whether 
it completely removes them. In any case the formalism for planar graph summation [7] is on 
firmer foundation, and is actually somewhat easier to implement, with a' > 0. 

The large A^ limit of gauge theory amounts to summing all the planar Feynman diagrams 
of perturbation theory. These diagrams are the a' — )■ limit of the planar open string multi- 
loop diagrams. Mandelstam [8] has given a remarkably simple and intuitive representation of 
the latter diagrams as path integrals over the worldsheets of the light-cone quantized string 
[9]. Each loop in a given multiloop diagram is represented as an internal worldsheet bound- 
ary whose beginning describes a breaking string and whose end describes two strings joining. 
For each fixed configuration of these internal boundaries, the path integral is a Gaussian in- 
tegral over the transverse string coordinates and the measure is precisely the natural lattice 
measure for the quadratic action. Summing over the number, location, and lengths of these 
internal boundaries with the same measure accomplishes the sum over planar open string 
diagrams [7]. 

While there is still some optimism that the large A^ limit of the maximal supersymmetric 

1 



(A/" = 4) gauge theory might be exactly solvable [3], it is doubtful that this will be possi- 
ble for the pure gauge theory underlying QCD. However, it might well turn out that the 
representation of the sum of planar diagrams by a lightcone worldsheet lattice system, as 
described in the previous paragraph, can be studied on a computer, perhaps using Monte 
Carlo algorithms. Such a numerical attack on QCD via its connection to open string the- 
ory could have strengths and weaknesses complementary to those of standard lattice gauge 
theory simulations. 

The basic framework for this approach to summing open string planar loop diagrams 
was set up by Giles and me (GT) over three decades ago [7] in the context of bosonic open 
string theory in 26 dimensional spacetime. The zero slope limit of these diagrams (modulo 
complications from the presence of the open string tachyon) would be 26 (not 4) dimensional 
gauge theory. Recently, I explained how to incorporate D3-branes into the GT formalism 
[10] in order to arrange zero slope limits that were four dimensional gauge theories, albeit 
coupled to 22 massless scalars corresponding to vibrations in the extra dimensions. There 
are mechanisms to decouple such massless scalars (e.g. see [11]), but the difficulties of the 
open bosonic string tachyon would remain. 

The open string Neveu-Schwarz model [12,13], restricted to even G-parity open string 
states (NS-I-), has no open string tachyon [14, 15], so the multiloop planar diagrams of the 
D = 4 version of NS+ would have a large A^ 4 dimensional gauge theory as a clean zero slope 
limit [16]. But since the standard lightcone quantization works only in the critical dimension 
[17-19], it is simpler to use the 10 dimensional version of the model, employing D3-branes to 
yield a four dimensional gauge theory. The 6 massless scalars can then be suppressed either 
by orbifold projections or by using nonabelian D3 branes [11]. 

The purpose of this article is to adapt the GT worldsheet lattice formalism to the Neveu- 
Schwarz model. This requires providing a viable lattice definition of the fermionic worldsheet 
field H^{a^ r) of that model. Each (transverse) component of if is a two dimensional Majo- 
rana spinor field. As always, putting fermion fields on a lattice involves difficulties. Besides 
the inevitable fermion doubling, for which there are a variety of remedies (Wilson fermions, 
staggered fermions, domain wall fermions), there is the difficulty that Grassmann path inte- 
grands do not have a probabilistic interpretation, a prerequisite for Monte Carlo methods. 
This last difficulty could be dealt with by integrating out the fermions, but that would 
sacrifice locality, rendering numerical simulations very costly. Fortunately, for the two di- 
mensional worldsheet, there is another option which we pursue here. This is based on the 
well-known fact that the physics of the critical two dimensional Ising model is that of a free 
Majorana Fermi field. The partition function of the Ising model is a sum over Ising spins 
Sij = ±1 with a positive definite Boltzmann factor, the exponent of which is local in the 
spins. By replacing each component of H with an Ising spin system, the lightcone lattice 
formalism for summing the planar diagrams of the Neveu-Schwarz model can be analyzed 
with Monte-Carlo methods. 

In this article we present and study a version of the Ising model which achieves this 
purpose. In Section 2 we recall the basic features of Onsager's solution of the model which 
employs a transfer matrix representation. We review the diagonalization of the bulk transfer 



matrix in terms of anticommuting spin matrices. In section 3 we turn to the issue of boundary 
conditions on an open strip, and show that a simple condition on the original Ising spins 
at the boundary produces the 1/2 integer modes required of the Neveu-Schwarz field H . 
Section 4 deals with the explicit construction of the eigenoperators of the transfer matrix. 
An important symmetry of the Neveu-Schwarz string is the so-called G-parity. In Section 
5 we identify the symmetry of the Ising model that becomes G-parity in the continuum 
limit. This is important for imposing the even G-parity restriction in lattice simulations. 
Unfortunately, imposing this restriction, in the most straightforward way, allows minus signs 
and/or nonlocality to creep back into the path summand, again posing potential difficulties 
for Monte-Carlo methods. In Section 6, we discuss the 3 open string vertex, which requires 
an operator insertion at the interaction point. It is argued that, provided the even G- 
parity restriction is maintained, these insertion factors can be taken into the exponent and 
interpreted as a modification of the worldsheet action. In Section 7 we give a brief discussion 
of the resulting representation of the sum over planar diagrams. Concluding discussion is in 
a final Section 8. 

2 Bulk properties of the Ising model 

We consider the two dimensional Ising model on an M x A^ lattice specified by the partition 
function 



z= Y, (^^^'^''''<+^^''''>^'^''' , (1) 



where 1 < i < M and I < j < N. Onsager's transfer matrix representation of Z is 



,\ MN/2 

Z = le-' -e-n {f\T{Mf\z) (2) 



J' ^ .~A^^^ !-( M\ - TTr.-^'/2 , ^x -J72 



(^■^e--fr(M) = ni^" 



a^e 




T{M) = exp^> afexp -> aM+i (3) 



tanh - = e"-^' , (4) 

where the states \i), \f) belong to an M-fold tensor product of 2-spinors, and they are 
determined by the boundary conditions at j = 1, A^ respectively. Here the o"f '^'^ are M 
independent sets of 2 x 2 Pauli spin matrices: 

{at,a':} = 2Sai„ [<,a|] = 0, for ^ ^ j . (5) 

An important property of the eigenvalue spectrum of the transfer matrix T{M) is an easy 
consequence of the representation (3): it is geometrically symmetric about 1. To see this 



note that if T has the eigenvalue T on the state \T), then it has the eigenvalue T^^ on the 
state e^^fe°'fc/^ nfc=odd'^fc 11/ "^f 1^)- ^^ shown by Onsager, the Ising model has a critical point 
when ^ = J. For the isotropic case J' = J, this occurs when sinh J = 1, or J = ln(l + ^/2). 
As usual, it is most convenient to replace the Pauli matrix dynamical variables with 
variables that anticommute for different i using the Jordan- Wigner trick 



rV'^ 



i-l 



fc=i 

{St, Sk} = SabSik , (7) 

whereupon the Onsager transfer matrix becomes: 

7-(jVf) = e-'^^^=^^'^^h-''^^^=^'^^^+^^' . (8) 

This choice for T imposes a particular set of boundary conditions at A; = 1, M appropriate for 
the Ising spins living on an open strip. If they lived on a cylinder, periodic or anti-periodic 
boundary conditions would be appropriate. We shall show that the boundary conditions 
chosen here lead to a half- integer moded fermion field in the continuum limit, which is what 
is needed to describe the Neveu-Schwarz model [12, 13]. 

It is straightforward to calculate the action of T by conjugation on the S variables: 

rS'^r-' = cjCi:Sl - icjs^Sl + isjc^Sl^^ - sjs^S'k_^, l<k<M (9) 

rSir-' = cjc^Sl + icjs^Sl-isjc^Sl^.-sjs^Sl^,, l<k<M (10) 

rs'^r-^ = c^s'^-is^sf (11) 

rsir-' = c^si + ts^si,, (12) 

where we have introduced the shorthand notation cj = cosh J, sj = sinh J, tj = tanh J, 
and similarly for J -^ C,- 

We see that T acts linearly on the S"s, which means that the problem of finding eigen- 
operators for T is one of linear algebra. The form of the right sides of (9), (10) suggest that 
an expansion of the S"s in eigenoperators will involve plane wave k dependence: 

S', = 5^A,e*'^ Sl = J2Bxe'^\ TiA,, B,)r^' = t,iA,, B,) . (13) 

A A 

We first consider the implications of the bulk recursion relations, focusing on a particular 
mode A. Suppressing the A subscripts on A, B, t, we find 

tA = cjc^A - icjs^B + isjc^Be^"-^ - sjs^Ae~''^ 
tB = cjc^B + icjs^A - isjCf-Ae^^^ - sjs^Be^^^ . 

The consistency of these two equations requires: 

\t - cjc^ + e-'^sjs^\^ = \cjs^ - e-'^sjc^\^ , (14) 



which simphfies to 

t^ - 2t{cjc^ - sjs^ cos A) + 1 = . (15) 

Since this equation is even under A — ?■ —A each of the two solutions, 



t± = cjc^ — sjs^ COS A =F \J {cjc^ — sjs^ cos A)^ — 1 , (16) 

belong to both left (A < 0) and right (A > 0) moving plane waves. This degeneracy is 
important in realizing definite boundary conditions. Also from the equation, it immediately 
follows that t_|_t_ = 1. 

The eigenoperators Aa-, with eigenvalue t^, connect eigenstates of T with different eigen- 
values. Because the state space is of finite dimension there is a maximum eigenvalue of T, 
Tmax, and a minimum eigenvalue: Tmin = T^l^ by the symmetry of the eigenvalue spectrum. 
Let \G) be the state with the maximum eigenvalue. Then we must have 

Aa\G) = 0, whenever t^ > 1 , (17) 



max • 



because if t^ > 1, then the eigenvalue of the state on the left would be taT^ax > T„ 
Assuming completeness of the eigenoperators, i.e. that monomials of eigenoperators acting 
on \G) span the whole state space, then the the state Ylt <i Aa\G) is the eigenstate with the 
minimum eigenvalue Tmm = T^ax nt„<i ta- It follows that T^^^x = Ilta>i ^«- 



Tma.= Y[ta. (18) 

Y ta>l 

Identifying In T^ax = —Eg the ground state energy, we can write this relation as 

ta>l ta>l 

where we have defined the energy created by an eigenoperator as a; = — In t. Then, since the 
eigenvalues t come in reciprocal pairs, the cj's come in opposite sign pairs w_ = — a;+, and 
we can write w± = ±u, with the convention that u > 0. 

In order that infinitely long waves (A — ?■ 0) cost zero energy, we see from (16) that the 
Ising system must be critical, that is ^ = J. In this critical case we can simplify 

t± = 



1 + (1-cosA)4t-^2(1- 


- cos A)sj + (1 — cos A)2sj 


(20) 


1 + 2s] sin^ - ^ 2 
It Asj . 


A 
s J sm - 


yl + ^sin'^ 


(21) 
(22) 



Once t is determined by consistency, we can then solve for B in terms of A: 

B = '-'i'^ + '~'''''^A^RA. (23) 



Notice that, according to the consistency condition, the modulus of the numerator of R is 
equal to the modulus of the denominator, so i? is a pure phase. Putting in the explicit forms 

for t±, we get 

-isjSi. sin A =F a/(cjC5 - sjs^ cos A)^ - 1 

ii± = TT ^T 7 • (^4) 

Specializing to the critical case, the form simplifies 



^^ ^ ^^x/2 -is J COS ( A/2) T sgn( A) yi + s] sin" (A/2) 

When we apply boundary conditions to an open chain in the next section, it will be necessary 
to consider linear combinations of left and right moving plane waves. Then we will take A, B 
as the coefficients of e^^^ and A' , B' as the coefficients of e~^^^ with A > in both cases. 
Then Eq.(24) stands and 

B' = R'A' (26) 

+isjs^ sin A T -J{cjc^ - sjs^ cos A)2 - 1 
""^ - W^^JJ^^^^) "~^"- ^'^ 

For the critical case these equations become 

^ ^ ^^x/2 -is J COS (A/2) T^/l + s] sin"(A/2) 

^, = ^-.V.-^^.cos(A/2)±Vl + .",sm-(A/2)^_^^^ ^^^ ^^^^ 

3 Boundary Conditions 

Now we turn to the issue of boundary conditions. The simplest open chain model, given by 
the transfer matrix of (8), will suffice for our purposes. The action of T on the first (S*") and 
last {S%i) spins (11), (12) must be made to agree with the bulk recursive formulas (9), (10) 
which were satisfied with our plane wave ansatz. This can be arranged by suitable definitions 
of 5*0 and 5*^^^, which do not appear in T. The bulk recursions read for A; = 1, M: 



y 

M+l ■ 



TSIT = cjCi-Sl - icjs^Sl + isjc^Sl - sjs^Sq 
TSIjT'^ = cjc^Sf^ + icjs^Slf - isjCi-Slj^^ - sjs^S] 

Agreement will occur if we impose 

{cj-l){c^St-zs^Sf) + sj{zc^S',-s^S',) = (30) 

{cj - l){c^Sl^ + is^Sl^) + sj{-ic^Sl^^^ - s^^lf+i) = . (31) 



Then plugging the plane wave ansatz 

S'k = Ae'^^ + A'e-^^^ SI = Be'^^ + B'e~'^^ (32) 

into these equations and rearranging leads to two different expressions for a defined by 
A' = aA: 

^ e'^(cj - l)(cg - JRs^) + sjjiRc^ - s^) 

e-i^(cj-l)ic^-iR's^) + sj{iR'c^-s^) ^ ^ 

^ ^2iMX (cj - l)(^cg + isg) + e'^sj{-ic^ - Rs^) 

{cj-l){R'c^ + is^)-e-^^sj{ic^ + R'si:)' ^ ' 

Remembering that R and R' = —R* have already been determined in (24), (27), we see that 
these two equations determine a and a quantization condition on A from the consistency of 
the two equations. We are particularly interested in the critical case ^ = J, for which the 
eigenvalue condition simplifies to 

g2iMA _ (e^^ - l)cj -e'^-l ^ _ l-icjtan{\/2) ^ ^^^^^^^ 

{e^ - l)cj + e^ + 1 l + icjtan(A/2) ' ^ ' 

We plot ^(A) in the critical case for three different values of J = .^ in Fig. 1. The graphical 
solution of (35) is shown in Fig. 2, for J = 1.5 and M = 5. This graph makes it clear that 
there are precisely M nontrivial eigenoperator solutions for any value of cj. The apparently 
(M+l)th solution A = vr is spurious. This is because, in this case, the left and right moving k 
dependence is identical (e*'^'^ = e"*'^'^) and the limit A — )■ vr implies C — ?■ —C and D' — )■ —D. 
Thus the corresponding eigenoperators, proportional to C + C" and D + D', vanish. For 
comparison, we show ^(A) and the solution of the eigenvalue equation for several noncritical 
values J < ^ in Fig. 3. 

Applying an eigenoperator to an eigenstate of T with eigenvalue Tq produces another 
eigenstate of T with eigenvalue tTo. We can say that the eigenoperator has increased the 
energy of the state by an amount u± = —lnt± = ±u. We show uj{X) for a critical case 
,^ = J in Fig. 4 and for a noncritical case in Fig. 5. The fact that the eigenvalue spectrum 
is non-degenerate guarantees that the set of eigenoperators is complete. 

The critical case allows a nontrivial continuum limit M — )■ oo: In that limit there are 
finite energy excitations which have A = 0(1/M). In this limit the eigenvalue condition 
simplifies even further to 

gi(2M+cj)A = _i + o{X^), for M ^ oo, AM fixed . (36) 

or A = {n+l/2)n/M, which is the mode quantization of the Neveu-Schwarz worldsheet field. 

4 Construction of Eigenoperators 

In this section we express the eigenoperators in terms of the S"s. We write the eigenoperator 

as 

M 

A = J2{SkUk + siVk), rAr-' = tA. (37) 

fe=i 
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Figure 1: 6{X) for the critical case ^ = 
(middle curve) and J = 1.5 (lowest curve) 



(highest curve), J 



Then we easily derive the recursion relations satisfied by the U, V: 

tUk = cjc^Uk- sjs^Uk+i + icjs^Vk-isjc^Vk^i, k = 2, 



tVk = cjc^Vk - sjs^Vk-i - icjs^Uk + isjc^Uk+i, 

tUi = C(Ui — SjS(U2 + icjs^Vi 

tVi = -isjJJi + cjc^Vi + isjC(U2 

tUM = CjC^Um - isjc^VM^i + is^VM 

tVu = c^Vm - icjs^UM - sjSjYm-x ■ 



k = 2,. 



M - 1 (38) 

M - 1 (39) 

(40) 
(41) 
(42) 
(43) 



The last four special cases will be included in the first two generic equations if we put 

isjVo = {cj - l)f/i, -isjUM+i = (cj - 1)Vm , (44) 

which specify the boundary conditions. As before, these relations can be solved with the 
plane wave ansatz: 



Uk = Ce'"^ + C'e 



k\ 



1/ „~ik\ 



ikX 



Vk = De'"^ + D'e 



\l „—ik\ 



(45) 
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Figure 2: Solution of (35) for J = 1.5. The straight hues are 2MA — 2n7r for 
n = 0,1, ... M — 1 for the case M = 5. Note that A = tt corresponding to 
n = M, though a solution of (35), is spurious as discussed in the text. 



'' ^ T ^ 



leading to 

D 

C 
D' _ D* 

together with the consistency condition 

t^ — 2t{cjc^ — sjs^ cos A) + 1 



iX, 







(46) 
(47) 

(48) 



which is identical to (15). Note that D/C is similar to B/A, but not identical to it. This is 
because the matrix we are diagonalizing is not hermitian. Of course, since the characteristic 
equation is the same for both eigenvalue problems, the eigenvalue spectrum will be identical. 
The boundary conditions read 

isj{D + D') = (cj-l)(Ce*^ + CV*^) 

which lead to 

C _ isj{D/C)-{cj-l)e^^ 
C isj{D*/C*) + {cj - l)e-^^ ' 

9 



(49) 
(50) 



(51) 




Figure 3: 9{\) for some noncritical cases: C, = 1.5 and J = 0.3,0.6,0.9,1.2 
from the highest to second lowest curves. The lowest curve is the critical 
case J = 1.5, included for comparison. The straight lines are 2MA — 2n7r for 
n = 0, 1, . . .M — 1 for the case M = 5. Their intersections with the ^(A) 
curves show how the solutions of the eigenvalue equation for the noncritical 
cases approach the solution for the critical case as J — )■ ^. 



and to the consistency condition 

,,j^, ^ jtsjiP/C) - jcj - l)e^')i-tsje-^' + jcj - 1){D*/C*)) 
{tsj{D*/C*) + {cj - l)e-^)((cj - 1){D/C) + isje^^) 
isjjl + \D/C\'^) + 2Re{e~'^D/C) + 2icjIm{e''^D/C) 
isj{l + \D/C\^) - 2Re{e-^^D/C) + 2icjlm{e-^^D/C) 
is J + Re{e~'^D/C) + icjlm{e~'^D/C) 
is J - Re{e-^^D/C) + icjlm{e-^^D/C) ' 



(52) 
(53) 
(54) 



The last line follows because one can show that C/D is a pure phase. Plugging the solution 
for t in the equation for D/C and a little manipulation shows that 



D 



-i\ 



—I 



± y^(cjc^"^^^js^7iosAp"^^ + isjs^ sin A 
cjS(^ cos A — sjc^ + icjs^ sin A 



(55) 
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Figure 4: The excitation energy uj{\) = — lnt+(A) for the critical case C, 
J = 1.5. 



±Jl + s'] sin2(A/2) + isj cos(A/2) 



cj(2 sin(A/2) + cos(A/2)) ' ^^^^ 

where the last line shows the simplifying critical limit .^ — ?> J. In this limit we also evaluate 

D' 



Re 



sj + cjim 



-i\^ 



-i\ 



C 

D 
C 



1 A 

cos — 

cj 2 



±W1 ' -2-"^ 



A 



St sm — \- Sj sm 



A 



^ / , /. 2 • 2 ^ 

sm — I ±\l I + s J sm — \- sj sm 



A 



Then the quantization condition on A simplifies dramatically to 

2iMX COS (A/2) — icj sin (A/2) 1 — icj tan (A/2) 

~ cos (A/2) + icj sin (A/2) ~ ~ 1 + icj tan (A/2) ' 

which is, of course, identical to (35). 



(57) 
(58) 



(59) 



5 G-Parity on the Lattice 

In the Neveu-Schwarz model the G-parity is a crucial symmetry concept. The fermionic 
raising and lowering operators b^ = b_r are moded with respect to i? = X]r>o^^-^^^ ^^ 
half odd integers. G-parity is defined by Ghj.G = —br, and G^ = 1. Defining fermion 
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Figure 5: The excitation energy w(A) 
J=l, C = 1.5. 



lnt+(A) for the noncritical case 



number as Np = J^ryo^-^K, G is essentially (— )^^. Because of the circumstance that the 
6's carry half odd mode number, we could also write (— )^^ = (— )^'^ which leads to technical 
simplifications in analyzing the model. The importance of the symmetry is that, since the 
interactions conserve G-parity, it is consistent to restrict the open string spectrum to be 
even under G-parity. Since the lowest odd G-parity state is a tachyon, and the lowest even 
G-parity state is a massless vector, the restriction to even G-parity removes the tachyon, 
leaving a massless sector that can describe a gauge boson. With this convention for even 
and odd, G = -(-)^^ = -(-)'^. 

There is no lattice analog of Np that commutes with the transfer matrix of the Ising 
model. However, there is a candidate for (— )^^, namely G = "/^Hfe'^fc' '^ith 77^ = 1. This 
operator, which anticommutes with Sf^, commutes with the transfer matrix. The Neveu- 
Schwarz model in D space-times dimensions requires D — 2 transverse components H and 
therefore D — 2 Ising systems, describable hj D — 2 sets of Pauli matrices cr^^'^ 1 ^ ~ 
1, . . . ,D ~2. Then the total G-parity operator will be 1] YlkA '^kA- It remains to fix the value 
oi 7] = ±1. According to the conventions of the Neveu-Schwarz model, we would like it to 
have the value —1 on the ground state of the system when it is critical. Since the ground 
state at general J has only been determined implicitly by rather complicated equations, it 
would seem a daunting task to evaluate the action of G on it. However, because G takes on 
only the values ±1, its value cannot change under the variation of a continuous parameter. 
Thus we can exploit simplifications in the state that occur when J — ?> at fixed ^ and M. 

As we can see from Fig. 3, the eigenvalue problem varies continuously as J decreases 
continuously from the J = .^ to J = 0. In the limit J — ?> the eigenvalue equation smoothly 
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approaches 



2iM\ „-2iA ^ ""-^ 



^,^M. ^ ^^.M ^^^__^ n=l,---,M. (60) 



Furthermore the eigenvalues t±{\) — ?■ e^'», independent of lambda, and the eigenoperators 
approach 



M , 

mrk 



A^^2^CJ2iSl±^Sl)sin——. (61) 



k=l 



M+1 



The A~ increase the eigenvalue of a state by the factor e^. Therefore in the limit the 
ground state of the system must be annihilated by all the A". Since the mode functions 
sin(n7rA;/(M + 1)) are complete, this limiting state must satisfy 

{S'k - iSl)\Go) = 0, forall A; = 1,---,M. (62) 

But af = iatal = 2iSlSl so it follows that 

al\Go) = 2iSl{-iS'k)\Go) = \Go), for all A; (63) 

HatlGo) = +|Go). (64) 

k 

By continuity we conclude that the ground state at all J < ^ will have this same eigenvalue. 
This conclusion depends on the equations determining the ground state not changing dis- 
continuously at some point. For example if the eigenvalue of an eigenoperator changed from 
t > 1 to t < 1 (from a; < to a; > 0) as J was varied, it would no longer annihilate the 
ground state and vice versa. But the character of the eigenvalue solutions doesn't allow this 
to happen, as is evident from Fig. 3. In particular note that the curve in Fig. 5 stays well 
away from u = throughout. We conclude that r] = —1, so the correct G-parity operator 
to use in describing the Neveu-Schwarz model is G = — YlkA ^kA- 

Finally, we work out the meaning of G-parity in the language of the original Ising partition 
function (1). To do this, focussing on a single Ising system, we consider the product of the 
appropriate factors of G with the transfer matrix (3): 



n ^ir = (e^' - e-^'r^l- n {'''''" + -y) -P U E -l-U^ ' (65) 

k k \ k ) 

compared to 

r = (e^' - e-^')-^/^ n {'■"" + -I'-'''") -P { ^ E -l^U^ I • (66) 

k V k ) 

Thus, inserting a G somewhere in the matrix element (2) has the effect on (1) of mul- 
tiplying Z by an overall factor (—1) and reversing the sign of one of the terms in the 
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sum over j, J' ^ji^j^A^kA'^kA )/2- Technically this is accomplished by inserting a factor 
—e^'^ ^kA^kA^kA J for a particular time slice /, in the summand of (1). To project onto even 
G-parity states, simply insert the factor (1 — e~"' ^ 
say it is to make the replacement 



kA^kA'^kA )/2. Or perhaps a better way to 



g^ J2kA '^kA'^kA /2 _^ sinh 




kA"kA 



(67) 



on at least one time slice of each open string propagator. This sinh factor is non local and 
also can be negative. The nonlocality would add to the computational cost of a Monte Carlo 
simulation. This is because the probabilistic criterion for retaining an update of the spin on 
a site /c, / requires knowledge of all the other sites on the time slice /, as well as those on 
neighboring time slices. However, the one-dimensionality of this nonlocality may help keep 
the cost manageable. The strict probabilistic interpretation of the summand of Z, which 
is the theoretical basis for Monte Carlo simulations, is marred by the fact that the sinh is 
negative for some spin configurations. However the spin configurations where the argument 
of the sinh is negative are strongly suppressed by the exponential factors on the many other 
time slices. One can then hope that those configurations will cause a minimal degradation 
of the convergence of the simulation. 

6 Open String Vertex 

A striking aspect of string theory is that interactions among strings are inherent in the 
nature of free string. This is because a single string can make a transition to two strings by 
simply breaking at a point. Using light-cone quantization of the free string [9] Mandelstam's 
interacting string formalism [8] provides the most concrete realization of this concept. As 
illustrated in Fig. 6, Mandelstam's three string vertex is simply a worldsheet path integral 
with the free string action, but for which the worldsheet fields live on a two dimensional 
domain corresponding to two strings joining ends to become a single string (or the time 
reversed process). Discretizing this domain on a worldsheet grid, Giles and I [7] clarified 
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Figure 6: Lightcone parameter domain for a three open string vertex. Light- 
cone time r is the horizontal axis and a is the vertical axis. 



the nature of the singularity induced by the string joining/splitting process for the bosonic 
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string in D space-time dimensions. Write T = aN, P^ + P^ = P+ = MaTo, so the 
diagram has dimensions P+ x T, and the associated lattice is M x A^. We found that, in 
the continuum hmit M, A^ — )■ oo, the 360° corner at the interaction point induced a behavior 
M~'^'^~^^/^^x Finite. For the critical dimension, D = 26, this scaling behavior accounts 
precisely for the {P^P^P+)-^/^ = {MaTo)-^/^{P^P^ /P+^y^/^ factor required by Poincare 
invariance. 

For the Neveu-Schwarz model, a fermionic worldsheet field, H^{z) = 'Yl,r^r^~^ with r 
ranging over half odd integers, is introduced. Of course, only the transverse components H^, 
k = 1, . . . ,D — 2 play a role on the lightcone worldsheet. The main thrust of this article 
is to represent these fields hj D — 2 independent critical Ising models on the worldsheet 
lattice. The contribution of the Ising degrees of freedom to the singular behavior at the 
interaction point can be inferred by realizing that a Majorana fermion is, roughly speaking, 
a half boson. More precisely, two Majorana fermion worldsheet fields can be interpreted, 
through bosonization, as a single bosonic worldsheet field. Thus the singular factor, including 
Ising and coordinate variables, should be M~^*^-^~^)/^^. If this were the end of the story, the 
Lorentz invariant critical dimension would be D = 16 + 2 = 18, not the Neveu-Schwarz 
critical dimension D = 10. 

The point is that the Neveu-Schwarz vertex is not simply the overlap represented by the 
diagram of Fig. 6, but there is also an operator insertion at the joining point [19]. In the 
continuum limit the insertion is just the density of the superconformal generators, H ■ x. 
Mandelstam showed that if this insertion is placed a distance e from the interaction point, 
then the amplitude ~ e^^'"^ as e — )> 0. More generally, if the operator insertion had conformal 
weight J, the singular behavior would be e~'^'^. 

We now translate these conclusions to a lattice calculation. Then there is no need to 
introduce e: one simply inserts the operator a lattice step or two away from the interaction 
site. Next, with lattice normalization the insertion operator itself would correspond to a"' 
times the continuum expression. For example, instead of x, one would insert xl —xl — )■ ax. 
Likewise, instead oiH which has delta function normalization, one would insert Sj — )■ a^^"^!!, 
which has Kronecker delta normalization. Thus putting e = a the continuum analysis of 
Mandelstam would lead to the behavior a"^ x a~"^/^ = a"^/^. Finally, in the lattice setup there 
is no reference to a, only to the number of lattice sites. Thus this estimate translates to the 
behavior 1/M"^/^ or 1/M^/^ for the case J = 3/2 of interest here. Putting this together with 
the result from the overlap leads to the overall scaling behavior 

/ -, N 3(D-2)/32+3/4 

^NS ~ [j^j (68) 

which is seen to give the correct Lorentz invariance power 3/2 for D = 8 + 2 = 10 in accord 
with the known properties of the model. This result has been obtained by translating the 
analysis done in the continuum theory into expectations for the results of a lattice calculation. 
It would be very nice to also see it from a direct lattice calculation, but we do not attempt 
that here. 

Now let us look more closely at the operator insertion from the Ising model point of 
view. We first need to decide which Ising model will describe the D — 2 = d fermi fields 
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^k'li ^ ~ ^)---)d- The most straightforward choice would be simply d decoupled Ising 
models as defined by (2), with each S^'^ built from the Pauli matrices following (6). Then 
to make the SkA for different A anticommute, one could define individual G-parity operators 
Ga = 11/0=1 "^fcA; ^^^ include a factor of n_B=i ^b in the definition of SkA'- 



A-l fc-l 



Sli = ^lill^^Il^iA, JWI. (69) 

B=l 1=1 

The trouble with this choice is that the insertion operator xaSa will have a residual nonlo- 
cality when expressed in terms of the original Ising variables. 

There is a better choice, which can be described as follows. Extend the index of S^'^ to 
the range k = 1, . . . , Md. Then identify the first d of these with Sia-, the next d with 5*2^, 
etc. Then relate the anticommuting S^.'^ to the Ising o"^'^ by the standard Jordan- Wigner 
transform (6): 

fc-i 
SV = (yV\{(^i, k = l,---,Md mil. (70) 

1=1 

The advantage of this version of the Jordan- Wigner trick is that the nonlocality of the 
insertion operator is subsumed in a factor which is proportional to the G-parity operator for 
one of the strings entering the vertex. When all strings entering the vertex are restricted to 
even G-parity, the nonlocality in the Ising variables disappears! Thus the insertion will be 
local in both the Pauli matrix and fermionic representations. This is very welcome, because 
it is the Pauli matrix form that will be more amenable to numerical analysis. On the 
unrestricted state space it is only the fermionic representation that is local. 

There is a price for this choice however. Since we demand that the transfer matrix 
expressed in terms of the fermionic spin variables be unchanged, the new Jordan- Wigner 
transform leads to a more complicated Ising model. To see how we plug the new Jordan- 
Wigner transform into the transfer matrix 

M c^i; Qz _i 7 V^d \^M-1 qz qy 

kA 



7~(Af) = e~*?EA=l Efc = l '^fcA'^feAg *'^EA=lEfc = l ^k + l,A^k. 
;c\^Md qy qz ,■ 7 Y^(*f-l)d QZ qy 



(71) 



The terms in the exponent of the first factor behave exactly as before: 

Md . Md Md 

k=l k=l k=l 

Md 

^^.r _^^J'^Md/2^-riinL\slsi ^ -Q(eJ' + ^Xe-J'). (72) 



fc=i 



16 



However the exponent of the second factor changes substantially: 



{M-l)d 



k=l 



J 



{M-l)d 



^■J Z^ ^k+d^k - t; / , 



z z 

"k+d^k 



k=l 



k+D-3 

n 

l=k+l 



a. 



Md 



(e^'_e-^')A^'^/2r(M) = W{e''/^ + ale-'' 



/2^ 



exp 



k=l 



J 



(M-l)d 



k+d-\ 



IT E ^^'^^i W^n (73) 



fe=i 



l=k+l 



Only for d = 1 does this reduce to the usual Ising model. The form of the Ising model 
corresponding to this transfer matrix is worked out in the Appendix. 

7 Summing Planar Loops 

Once the three open string vertex has been established as in the previous section, the com- 
plete perturbation series is determined [8, 19]. A generic planar multiloop diagram in the 
series is the worldsheet path integral using the free open string worldsheet action, but with 
the worldsheet variables living on a domain such as depicted in Fig. 7. 



p+ 



Figure 7: A planar multiloop lightcone interacting string diagram. The hori- 
zontal lines form the boundaries of the propagator worldsheets for intermediate 
open strings. This diagram has 7 loops and 5 external strings. 



In order to digitize the sum over planar diagrams, we fix the overall dimensions of the 
domain to P^ x T, and then set up a worldsheet grid of dimensions M x N, with T = aN 
and P~^ = oTqM [7]. The discretized worldsheet action is constructed so that the internal 
horizontal lines shown in Fig. 7 represent open string boundaries. For the Neveu-Schwarz 
model the worldsheet variables are the d = D — 2 transverse coordinates x{(t, t) — )■ xj and 
the d fermionic fields y/aH{a, r) — )■ Si where the S^ are taken to be the Jordan- Wigner 
transformed Pauli spin matrices of d independent Ising models. For each end of a horizontal 
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line, which depicts the breaking or joining of open strings, there is a factor of the open 
string couphng g and also the operator insertion 5/ ■ (xj'^ —^l) as explained in the previous 
section. The precise location of this insertion is somewhat flexible, as long as it is within 
one or two lattice steps from the end of the horizontal line. For definiteness, we will always 
place it on the longest string participating in the vertex, with i the spatial location of the 
horizontal line, and j — 1 or j ' + 2 the time of the end of the horizontal line, with the choice 
determined so that the insertion lies completely on the longest string participating in the 
vertex. 

Next we turn to G-parity restrictions. The vertex obtained in the previous section, and 
adopted in this section, conserves G-parity: it connects 3 even G-parity open strings with 
each other or 2 odd G-parity strings to an even one. We would like to restrict the open 
string states to even G-parity only. When a diagram involves one or more loops, it is not 
sufficient to restrict the external states to even G-parity, because a pair of odd G-Parity 
states can be produced from an even G-parity state. Thus each internal propagator in a 
multiloop diagram such as depicted in Fig. 7 must include a projector {1 + G)/2 onto the 
even G-parity sector. 

As we have discussed in the previous section, the presence of even G-parity projectors 
throughout the diagram introduces nonlocality into the worldsheet dynamics. The projectors 
also produce negative contributions to the path integrand. These could potentially lead to 
inefficiencies and inaccuracies in the applications of Monte Carlo algorithms to this system. 
However, there are beneficial aspects of the presence of the projectors. One, already men- 
tioned, is that restriction to the even G-parity states renders the nonlocality in the relation 
between the S"s and the a's harmless. Thus the operator insertions, necessary to describe 
the Neveu-Schwarz model, will have a local representation in terms of the original Ising vari- 
ables. Even so, the insertions are awkward because they are indefinite in sign and have no 
obvious interpretation as terms in the worldsheet action. 

It would be nice if these factors could be taken into the exponent where they would 
become a modification of the action^. If we do that and then expand the exponential, 
the effect would be to replace the desired insertion with a sum of all possible powers of 
the insertion, including a term with no insertion at all. The higher powers are innocuous 
because they will either renormalize the zeroth and first powers or introduce operators of 
higher conformal weight which would be suppressed relative to the zeroth and first powers 
in the continuum limit. It is the zeroth term that poses the difficulty. In the continuum 
limit it would dominate over the desired linear term. On the other hand it does not couple 
3 even G-parity open strings together. Thus, if the even G-parity restriction is enforced 
throughout, this troublesome term will not contribute. So by including the projectors, we 
enable the interpretation of the operator insertions as modifications of the action. 



^Whcn H is represented as a Grassmann variable this procedure would have the bizarre effect of adding 
Grassmann odd terms to the action. However in the Ising model representation it would simply mean adding 
terms linear in the spin to the action-no weirder than an external magnetic field! 



8 Concluding Remarks 

In this article we have extended the GT lattice lightcone string formalism to include the 
Neveu-Schwarz model, representing the Fermi fields via Ising spin systems. That represen- 
tation has the virtue of avoiding the minus signs inherent in a description of fermions as 
Grassmann fields. If there were no need to include operator insertions at the interaction 
points and no need to make the even G-parity restriction, the formalism would be ready for 
immediate analysis via Monte Carlo simulations. 

Alas, that would give neither the Neveu-Schwarz (NS) model nor its even G-parity pro- 
jected descendent (NS+). Indeed, the simulations would be dominated by tachyon effects, 
most likely yielding the same sort of almost trivial outcome as the bosonic string [20], namely 
the sum of planar open string diagrams would produce the propagator of a free closed string. 
What makes the NS+ model dynamically interesting is precisely what makes its Monte Carlo 
analysis problematic: the operator insertions in the NS model are of indefinite sign and the 
projectors necessary to describe the NS+ model are nonlocal. Intriguingly, for the NS+ 
model, we argued that the restriction to even G-parity provides a resolution to the operator 
insertion difficulty. In its present status, all of the sign and nonlocality problems of the for- 
malism reside entirely in the projection procedure. We think it is some progress to be able 
to attribute all of the nonlocality and negativity to such a well-defined source, in which they 
may turn well out to be relatively benign. In any case, we suspect that a better formulations 
will be found to circumvent any difficulties that remain. 

Turning to more mundane matters, we recall that our main motivation for keeping a' > 
was to mitigate the need for uncontrolled counterterms to maintain Lorentz invariance in 
the loop expansion. While this should effectively deal with the usual field theoretic UV 
divergences, it does not by itself take care of the potential need for the worldsheet contact 
counterterms studied, for example, in [21]. This issue clearly needs further study. However it 
is comforting to know that the fundamentally stringy regularization identified in [22] (GNS), 
seems to be at least partially realized by the GT lightcone lattice. The GNS mechanism, 
works by regarding each open string loop on the worldsheet as an emission or absorption 
of a closed string in the vacuum: giving that closed string a nonzero momentum regulates 
the worldsheet divergence. By studying the one loop open string self energy, we have found 
that discretizing P"*" = Ma, which is the discretization of GT lightcone lattice, has a similar 
regulating effect as p 7^ 0. For this example it guarantees that the open string gluon remains 
massless! If this continues to happen in more complicated multi-loop processes, we will be 
in business. 
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A Alternate Description of d Ising Spin Systems 

Here we work out the Ising model corresponding to the transfer matrix (73). The way to do 
this is to evaluate the transfer matrix elements between basis states which are eigenstates 
of the o"^. Then T^ can be expressed as sums over these eigenvalues. However it is more 
convenient to factor T into d = D — 2 factors, and work out the matrix elements of each 
factor. So write 



ieJ' - e-^V^^/2 



nM) = n n(^''''+-; 



M 



'm+(k-l)d^ 



.-J'/^] 



m=l '- k=l 
M-l 



k=l 



J 



Y\ exp < — (Tl , i,^a, 



2~ m+kd'-'m+{k-l)d J_ J_ '^l 

l=m.+l+(k-l)d 



m—l+kd 

n 



(74) 



and we evaluate the matrix elements of each T^, where 



T 



,J'/2 , ^x 



+ (t: 



m+(fc-l)d^ 



M 

n( 

k=l 

M 

k=l 



-J'/2^ 



M-l ^ J m-l+kd ■. 

n ^^P ] 2^m+kd(^m+{k-l)d H ^^ \ 

fc=l *- l=m+l+(k-l)d ^ 



"' - <+(.-l).e-^'/2, 



M-l 



VI— L y- 

n 

k=l ^ 



cosh£^..W™..-.. 



m—l+kd 



TT af sinh 



J 



l=m+l+{k-l)d 



(75) 



Notice that in the last expression no a^ appears more than linearly. Thus we can use the 
identities 



(aVk) 



I -a' a 



{aVW) 



l + a'a 



2 ' V . . / 2 

where \a) is an eigenstate of a^ with eigenvalue a. Then a few minutes thought leads to 

M 

l)d^m+{k-l)d 



(76) 



Cjt M 

^ k=l 



M-l 

n 

fc=i 



m—l+kd _, 

n '-^ 

l=m+l+{k-l)d 


— cosh — 

2 




m—l+kd 


+0'm+kdO'm+{k-l)d 
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I 


=m+l+{k-l)d 



i^^sinh:^' 

2 2 



(77) 



The matrix element of the full transfer matrix is then 

d 

(Ml n ^™l W) = E (Ml^il{^i})(W}l^2|{^2}) ■ ■ ■ (M-3}|T.| W) . {n 



m=l 



{R}} 
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In effect we can thinlc of each actual time slice as d coincident time slices, so that the 
computational scale of this new version of the Ising model would correspond to an M x dN 
lattice as compared to an M x A^ grid of the simple Ising model. And of course as complicated 
as this formulation seems, it has identical physics to d decoupled Ising models. The only 
reason for suffering these extra complications is that in the application to the interacting 
NS+ model, the insertion operators will be local in the new representation, whereas they 
would be nonlocal in the original formulation. 

References 

[I] A. Neveu and J. Scherk, Nucl. Phys. B 36 (1972) 155. 

[2] J. Dai, R. G. Leigh and J. Polchinski, Mod. Phys. Lett. A 4 (1989) 2073. 

[3] J. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), hep-th/9711200. 

[4] G. 't Hooft, Nucl. Phys. B72 (1974) 461. 

[5] K. Bardakci and C. B. Thorn, Nucl. Phys. B 626 (2002) 287 [arXiv:hep-th/0110301]; 
C. B. Thorn, Nucl. Phys. B 637 (2002) 272 [Erratum-ibid. B 648 (2003) 457] [arXiv:hep- 
th/0203167]; S. Gudmundsson, C. B. Thorn and T. A. Tran, Nucl. Phys. B 649 (2003) 
3 [arXiv:hep-th/0209102]. 

[6] D. Chakrabarti, J. Qiu and C. B. Thorn, Phys. Rev. D 72 (2005) 065022 [arXiv:hep- 
th/0507280]; D. Chakrabarti, J. Qiu and C. B. Thorn, Phys. Rev. D 74 (2006) 045018 
[Erratum-ibid. D 76 (2007) 089901] [arXiv:hep-th/0602026]. 

[7] R. Giles and C. B. Thorn, Phys. Rev. D 16 (1977) 366. 

[8] S. Mandelstam, Nucl. Phys. B 64 (1973) 205; 

[9] P. Goddard, J. Goldstone, C. Rebbi and C. B. Thorn, Nucl. Phys. B 56 (1973) 109. 

[10] C. B. Thorn, Phys. Rev. D 80 (2009) 086010 [arXiv:0906.3742 [hep-th]]. 

[II] C. B. Thorn, Phys. Rev. D 78 (2008) 106008 [arXiv: 0809. 1085 [hep-th]]. 
[12] A. Neveu and J. H. Schwarz, Nucl. Phys. B 31 (1971) 86. 

[13] A. Neveu, J. H. Schwarz and C. B. Thorn, Phys. Lett. B 35 (1971) 529. 

[14] S. Mandelstam, private communication, 1971. 

[15] F. Gliozzi, J. Scherk and D. I. Olive, Phys. Lett. B 65, 282 (1976); Nucl. Phys. B 122 
(1977) 253. 

[16] C. B. Thorn, Phys. Rev. D 78 (2008) 085022 [arXiv:0808.0458 [hep-th]]. 

21 



[17] P. Goddard and C. B. Thorn, Phys. Lett. B 40 (1972) 235. 

[18] P. Goddard, C. Rebbi and C. B. Thorn, Nuovo Cim. A 12 (1972) 425. 

[19] S. Mandelstam, Nucl. Phys. B 69 (1974) 77. 

[20] P. Orland, Nucl. Phys. B 278 (1986) 790. 

[21] M. B. Green and N. Seiberg, Nucl. Phys. B 299 (1988) 559. 

[22] P. Goddard, Nuovo Cim. A 4 (1971) 349; A. Neveu and J. Scherk, Nucl. Phys. B 36 
(1972) 317. 



22 



